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Hankel hyperdeterminants and Selberg integrals 
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2, France 

Abstract. We investigate the simplest class of hyperdeterminants defined by Cayley 
in the case of Hankel hypermatrices (tensors of the form Ai^i^ = /(ii + b + ’ ■ • + *fe))- 
It is found that many classical properties of Hankel determinants can be generalized, 
and a connection with Selberg type integrals is established. In particular, Selberg’s 
original formula amounts to the evaluation of all Hankel hyperdeterminants built 
from the moments of the Jacobi polynomials. Many higher-dimensional analogues of 
classical Hankel determinants are evaluated in closed form. The Toeplitz case is also 
briefly discussed. In physical terms, both cases are related to the partition functions 
of one-dimensional Coulomb systems with logarithmic potential. 


1. Introduction 


Although determinants have been in use since the mid-eighteenth century, it took almost 
one hundred years before the modern notation as square arrays was introduced by Cayley 
1^. Then, it was not long before Cayley raised the question of extending the notion 
of determinant to higher-dimensional arrays (e.g., cubic matrices and proposed 

several answers, under the name hyperdeterminants 1,0- 

The most sophisticated notion of hyperdeterminant has been the object of recent 
summarized in the book []T|. However, the simplest possible 


investigations, 

generalization of the determinant, defined for a fctli order tensor on an n-dimensional 
space by the fc-tuple alternating sum (which vanishes for odd k) 

1 


Detfc(A) = — ^ e(ai)---e(crfc)||A^pi)...„,(i) (1) 

Cl.-'-jO-feeSn i=l 

has been almost forgotten. After the book by Sokolov which contains an exhaustive 
list of references up to 1960, we have found only [^, [T^, |^. These references 


contain evaluations of a few higher dimensional analogues of some classical determinants 
(Vandermonde, Smith, ...). However, the analogues of Hankel determinants do not 
seem to have been investigated. 

In this article, we shall compute the hyperdeterminantal analogues of various classical 
Hankel determinants. The elements of these determinants will in general be 
combinatorial numbers or orthogonal polynomials. Our main technique will be the 
use of integral representations. We shall see that the relevant tool is Selberg’s integral 
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and its generalizations, mainly in the form given by Kaneko. More that this, we can 
say that the knowledge embodied in Selberg’s formula and its limiting cases amounts 
to a closed form evaluation of all Hankel hyperdeterminants built from the moment 
sequences associated with the classical orthogonal polynomials. 

A more general class of hyperdeterminants is given by the partition functions Zn{f3) 
of log-potential Coulomb systems, when (3 is an even integer. When the particle- 
background interaction does not lead to a Selberg integral, the partition function can 
usually be evaluated in a more or less closed form only for (3 = 1,2,4 (as a Pfaffian or 
a determinant). The evaluation of general hyperdeterminants is of course much more 
difficult, but their simple transformation properties leave some hope that at least some 
of them may be evaluated by higher-dimensional analogues of the algebraic techniques 
working for Hankel determinants. 

The consideration of the partition functions of similar Coulomb systems with the 
particles conhned on a circle suggests immediately the following dehnition. A tensor 
called a Toeplitz tensor if 

= /(h + ■ • • + 4 - ji-jfc) • (2) 

Indeed, when (3 is an even integer, the partition function turns out to be a Toeplitz 
hyperdeterminant. 


2. Hankel hyper determinants 


Let be a tensor of order k and dimension n. The tensor A is said 

to be a Hankel tensor if = /(A + ■ ■ ■ + ik)■ 

Let us now fix some sequence c = (cn)n>o, and consider the hyperdeterminants 

^n\c) = Det2fc(Cj^+...+j2fc)o<ip<n-l ( 3 ) 


as dehned by formula (|l]), in which G„ is the symmetric group and e(cr) the signature 
of a permutation a. For k = 1, this is an ordinary Hankel determinant. For n = 2, it is 
easy to derive the expression 

ofM = lf;(-i)‘(^Ac.Q.-i (4) 

i=o ^ ' 


whose right-hand side is well-known in classical invariant theory (it is one-half of the 
apolar covariant of the binary form f{x,y) = with itself, see [pT| ). 

The case c„ = n! will be used as a running exemple throughout this paper. Using (^, 
we can give our hrst illustration of a higher-order determinant 

2k 


D 


(k), 


1 


c = 


5 ^(- 1 )‘( 2 «:)! =-( 2 *:)! 


( 8 ) 


i=0 

which will provide a check for the general case. 

Let now /i be the linear functional on the space of polynomials in one variable 
such that p{x"') = Cn- We extend it to polynomials in several variables by setting 
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rpTTl-ri I — 

) ~ ^rni 

determinant 


Cn 


^j) = 


Then, using the expansion of the Vandermonde 
^ e(cT)cT«-^a;”:^■•a:2) (6) 


i>j cr^Gn 

where a permutation a acts on a monomial by sending each Xi on a;o-(j), it is easily seen 
that 


(C) = - 7 /‘«(A“(^)) 


n\ 


(7) 


Expanding each factor {xi — Xjf^ by the binomial theorem, we obtain 


2k 


m 




n^apCM) 

p=l 


( 8 ) 


M=(mij) i>j 

where M runs over all strictly lower triangular integer matrices such that 0 < mij < 2k, 
\M\ = and 

n p—1 

ap{M) = 2k{p - 1) + ^ ruip-'^ rupj . (9) 

2=p+l j=l 

This extends and provides a faster algorithm than the definition. 

Now, if /i is a measure on the real line, then 


= — A‘^^{x)dp{xi) ■ ■■dii{xn ). 


n\ 


( 10 ) 


When k = 1, this is a well-known formula due to Heine [Q. For arbitrary k, the integral 
can be evaluated in closed form in many interesting cases by means of Selberg’s integral 
formula which gives, for 


Sn{a, h, k) = 


\A{x)\^'^\{xt\l-x,f-^dx, 


2=1 


the value 


n—1 


Sn{a,h, fc) = JJ 


T{a + jk)T{b + jk)Ti{j + l)k + l) 
^{ 0 - +b + {n + j - l)k)T{k + 1) 


( 11 ) 


( 12 ) 


The formula is valid, when defined, for complex values of k as well, but it is interesting to 
observe that all its known proofs start with the assumption that /c is a positive integer, 
and then extend the result by means of Carlson’s theorem (see, e.g., H). 

Hence, Selberg’s integral computes precisely the Hankel hyperdeterminants D^\c) 
when the Cn are the moments of the measure dp{x) = l[o,i]a;““^(l — x)^~^dx (the Beta 
distribution). The orthogonal polynomials for this measure are, up to a simple change 
of variables, the Jacobi polynomials — 2x). Hence, we can as well compute 

Dn\c) for the moments of the ordinary Jacobi polynomials, and their limiting cases 
(Laguerre and Hermite). The appropriate variants of Selberg’s formula are listed in |3(] 
17.6. 
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When the measure dfi{x) does not lead to a known variant of Selberg’s formula (such 
as Aomoto’s and Kaneko’s generalizations, which are discussed below), it is sometimes 
possible to evaluate the hyperdeterminant of order 4 {k = 2) from the knowledge of the 
scalar products (P„, P^) of the corresponding (monic) orthogonal polynomials and their 
derivatives. Indeed, it is classical (see |^) that 

A(a:)^ = det{Pi-i{xj)\P'_^{xj)) (13) 

where, in the right-hand side, we mean the 2n x 2n-matrix with i = 1,... ,2n, 
j = 1,... ,n, whose first n columns are the Pi and the last n ones their derivatives. 


Using one of de Bruijn’s formulae (see A. 18.7), we can write the hyperdeterminant 

as a Pfaffian 


Det4(ci+j+fc+z)|o = pf(My)| 


|2n-l 

0 


where M is the skew symmetric matrix such that 


(14) 


(15) 


M,, = {p„p;)= / p,{x)p;{x)dp{x) 

J a 

a i < j. For example, if dfi{x) = e~^dx is the Laguerre measure on [0, cx)), whose monic 
orthogonal polynomials are given by the generating series 

-t)” 




[x) 


n>0 

an easy calculation gives 


n\ 


xt 

et-i 

1-f 


f _ 1 'vm+n 


rs"" = 


n,m>0 


nimi 


(l-s)(l-st) 


which leads to 


Det 4 ((i + j + k + 


(16) 


( 17 ) 


(18) 


with = 1 if z > j and 0 otherwise. The calculation of the Pfaffian is straightforward, 
and we obtain 

2n-l 

Det4((z-|-J-|-fc-|-/)!)I q (19) 


i=0 


a special case of the general formula (P^D derived below from Selberg’s integral. 

For later reference, let us recall that the Hankel determinants Pn\c) are the products 
of the squared norms of the monic orthogonal polynomials P„, and that, more generally, 
the shifted Hankel determinants associated to the shifted sequences 

(r) . 1 

Cn = Cn+r cLie giveu by 


B<‘>(c) = det({Vf;,P,))|y‘, (20) 

Similarly, the shifted hyperdeterminants of order 4 can be reduced to Pfaffians 

OLfM = pf ((-!)*«> (VP., . (21) 
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Finally, let us remark that when the measure can be written in the form dp{x) = 
Ce~^^^^^dx, the integral (p!0D represents the partition function of a one-component log- 


potential Coulomb system on the line, evaluated aX (3 = 2k (see, e.g., [j^). It is a 
common feature of most of these systems that the partition function may be evaluated 
in closed form only for fd = 1,2,4 (the case 13 = 1 can be formulated in terms of Pfaffians 
of bimoments of skew-orthogonal polynomials, see |^, and will not concern us here). 
Similarly, the partition function of a one-component Coulomb system of n identical 
particles conhned on the unit circle has the general form 

N 


Zn{f3) = 




.-dvizj) 


( 22 ) 


For (3 = 2/c, this is, up to a scalar factor, the hyperdeterminant of the Toeplitz tensor 
associated to the bi-in£nite sequence 


Cr?, — 


1 

2'xi 


^-n^-yviz) 


dz 


(23) 


l«l=i 


As above, the knowledge of the appropriate orthogonal polynomials allows one to 
evaluate the determinant {k = 1) and sometimes the 4-fold hyperdeterminant. 


3. Examples involving combinatorial nnmbers 

The evaluation of Hankel determinants built on classical sequences of combinatorial 
numbers arises in many contexts (see p0| , ^ and references therein). Also, recent work 
on the theory of coherent states has led to the discovery of integral representations 
of many such sequences, in a form directly relevant to the evaluation of their Hankel 
hyperdeterminants 1^^ ^ . 


3.1. Factorials and T-functions 

As a warm-up, let us start with the already considered sequence 

poo 

Cn = n\ = F(u + 1) = x‘^e~^dx. 

Jo 


Then, using the Laguerre-Selberg integral (see |^, (17.6.5)) 

r(l + 7 + J7)r(a +jq) 


LSn{a,'y)= / |A(a;)p'^ ^^dxj = Y\_ 

J{0,oo)" 

we hnd 

n—1 


r(i + 7) 


F(u, k) := Y[{k + kj)\{kj)\ 


j=0 


thus recovering the classical evaluation Dn\c) = [l!2! ■ ■ ■ (n — 1)!]^ (see, e.g., 
the Hankel determinant. At no more cost, we can take 

Cn = F(n -|- 1 -|- a) 


(24) 

(25) 

(26) 

) of 


(27) 
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corresponding to the measure dfi{x) = x°‘e~^dx and obtain 

-l n—1 

^ n (*; + + l+jk) (28) 

j=0 

(k) 

which, for a = r a. positive integer, gives the shifted hyperdeterminant Dn-,r for factorials. 


3.2. Catalan numbers 


Here, we take Cn = Cn = C^n)- Penson and Sixdeniers give in an integral 

representation of these numbers 

1 ri 


C„,= 


27r 


X 


X 


—dx = -B(n + i, |) 


TT 


(29) 


The Hankel hyperdeterminant can be written in the form 

o2fcn(n—l)+n 

■Snd.I.fc) 


n'.TT" 


2 fcn(n-i)-n ^ kj)\{2kj + l)!!(2A;j - 1)!! 


n\kV 


n 

i=o 


(1 + k{n + j - 1))! 


(30) 


where (2n + 1)!! = 1 ■ 3 ■ ■ ■ (2n + 1). The Hankel hyperdeterminants of shifted Catalan 
numbers can be obtained similarly. Indeed, replacing c„ by Cr+n leads to another Selberg 
integral 


(y2kn{n—l)-\-n{2r-\-l) 

= -773^-S«(>- + 5,1.*;) 


mvr" 


(31) 


In the case fc = 1, the Hankel determinants of shifted Catalan numbers have been 
computed by Desainte-Catherine and Viennot [Q with the aim to enumerate the Young 
tableaux whose columns consists of an even number of elements and have height at 
most 2n. A natural question is to hnd what is counted by hyperdeterminants of shifted 
Catalan numbers. 


3.3. Central binomial coefficients 


In [321, Penson and Solomon give the representation 


2n 

n 


= — X 
TT ./n 


, 4'" 

*[a;(4 — x)]~^^'^dx = — B(?7, + i, D 

TT 


SO that, for c„ = (^"), 


^/cn(n—l)+nr 


D'nHc) = ' SJlkk). 


n'.TT" 


Similarly, the shifted hyperdeterminants are given by 

Akn{n—l)^nr 

Di^l = -r-- Snir + I, i, k) . 


n'.TT'' 


(32) 


(33) 


( 34 ) 
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Here, we find in that 

= J_ r x^e-^/^x-^/^dx . (35) 

Setting X = Ay and using the Laguerre-Selberg integral (p5|), we obtain the shifted 
hyperdeterminants as 

Di^c) = k) 

p§fcn,(n—l)+rn 

=--jTp^- YliH^ + j))K‘^ikj +r)- 1)!!. (36) 

j=0 


3.5. Bell numbers and polynomials 


We now take Cn = bn{a), the (one-variable) Bell polynomials, dehned by 

n 

6o(a) = 1 and bn{a) = S'(n, k)a^ (37) 

k=l 

where the S{n, k) are the Stirling numbers of the second kind (so that &n(l) are the Bell 
numbers). These are the moments of the discrete measure 

dpa{x) = e~''^^6{x - k) (38) 

A:>0 

for which the Charlier polynomials are the orthogonal system (cf. [|^). The monic 
Charlier polynomials Cn\x) satisfy 

(Ct‘>,Cy)=a"n! (39) 


whence the classical evaluation of the Hankel determinants 

n—1 

p){l) ^ _ 

j=0 


(40) 


However, no analogue of Selberg’s integral is known for the measure dpa- So, the best 
that we can do is to evaluate the fourth-order hyperdeterminants by means of formula 
(|T^). To this aim, we need the scalar products Cm'), which can be easily obtained 
from the generating function 

7/” 

C{n, x; o) = 5 ^ = e-“(l + «)“ ■ (41) 

n>0 


Taking the scalar product of this expression with that 


(cy,c£>') = 


(- 1 ) 


m—n+l g^m! 

m—n 

0 


dx 

ii m> n 
otherwise. 


(42) 
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It would remain to find a closed expression for the Pfaffian (|^. The hrst values are 


D^\c) = a(l + 6a) 

L)f)(c) = 8a%l + 24a + dSa^ + 90a^) 

Df\c) = 1728a®(l + 60a + 360a2 + 2080a^ + 2415a^ + 2100a® + 2100a®) 


It appears that Dn'^ (c) is always divisible by Dn'^ (c) (which is not true for general Hankel 
hyperdeterminants). For k > 2, we can compute the hrst polynomials: 

D^ 2 \c) = a(l + 30a + 60a^) 

Df\c) = 32a®(l + 240a + 3285a2 + lOOSOa^ + 61425a^ 

+ 56700a® + 37800a®) 

D^^\c) = a(l + 126a + 840a2 + 840a®) 

Df{c) = 128a®(l + 2184a + 134505a2 + 1952370a® + 22027950a^ 

+ 99542520a® + 189552825a® + 246673350a^ 

+ 130977000a® + 43659000a^) 

It is also interesting to observe that the shifted determinants Dn)r{c) can be expressed 
as Wronskians 

Dd)(c) = hr+u • • •, Wi)(«) • (43) 

This identity follows immediately from the recursion 

hn+i{a) = a[hn{a) + &'„(a)] (44) 

and is not of the same nature as the Karlin-Szego-type identities like (^^ below, 
discussed in section in which the shifted Hankel determinant of order n is expressed 
in terms of a Wronskian of order r. Here, also, Zli^r(c) is always divisible by Dn\c), 
as can be checked from Table 1. In this case, the explanation is simple: it follows from 
(pop that 

£>«(c) = det((VP,, P;))D<i\c) (45) 

where P* is the adjoint basis of P,. The ratio Pn:r(c)/Dlt\c) is therefore the determinant 

^ (r) 

of the operator Xn ^ = Mr o vr^ where Mr is multiplication by x'^ and 7r„ the orthogonal 
projection on the subspace spanned by Pq, ..., Pn-i- The matrix elements of Xn'^ can 
be read directly on the three-term recurrence relation of the monic polynomials, by 
iterating it, if necessary, to express as a linear combination of the Pj’s. The matrix 
element {x^Pi, P*) is then equal to the coefficient of Pj in this expression. 

For example, if the Pn = Ci^\x) are the monic Charlier polynomials, the three-term 
recurrence is 

xPi = Pj+i + {i + a)Pi + iaPj_i (46) 

so that 

x^Pj = Pj+2 + (2i-l-l-|-2a)Pj+i-|-(a^-|-a-|-4ai-l-i^)Pi-l-ia(2i — \p2a)Pi-i-\-i(i — l)aPi_2(47) 
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Table 1. The first values of 


n\r 

0 

1 


2 




1 

1 

a 


a 

f a2 



2 

a 

a? 


a® 

(2 2ci 

a^) 


3 

2a^ 

2a® 


2a® (6 + 6a + 3a^ + a®) 


4 

12o® 12a 

10 

12a®0(24 + 24a + 12a2 

+ 4a® + a^) 

Table 2. 

The triangle T^j 



k\j 

1 

2 

3 


4 

5 

6 

1 

1 







2 

1 

6 






3 

1 

30 

60 





4 

1 

126 

840 


840 



5 

1 

510 

8820 

25200 

150120 


6 

1 

2046 

84480 

526680 

831600 

332640 


and for r 


2 and n 



3, the matrix is 
’ a + a + 2a^ 

1 “t“ 2(2 1 5(2 “h (2^ 

1 3 + 2a 


2a^ 

6a + 4a^ 

4 9(2 -|- (2^ 


(48) 


whose determinant is 6a^ + 6a^ + 3a® + a®. This is the valne at x = 0 of the 2-Wronskian 
W{C^\C^^^){6) of Charlier polynomials. There is a general formnla (apparently new) 




4 “) 


l!2!---(r- 1)! 


(49) 


which will be derived in Section 

The seqnence D^\c) of bidimensional hyperdeterminants gives rise to an interesting 
triangle of integers T^j dehned by 

k 

of(c) = 5^r«o^ (50) 

i=i 

whose hrst valnes are given in Table 2. The main diagonal is given by {2k + l)!/fc!, and 
the second colnmn is 2^^“^ — 2. It is not difhcnlt to give a generating fnnction for these 
nnmbers. According to (^), if we know the exponential generating fnnction 

g{x) = J2cn^ (51) 

n>0 


then 


1 

2 




k>l 


x2k 

M! 


1 

2 


g(x)g(-x). 


(52) 
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1 ^2fc 1 

- + ^ dP(c)—, = -exp|o(e" + e'" - 2)]. 


fe>l 


(2k)\ 


10 


( 63 ) 


Since the hn{a) are the moments of a discrete measure, according to the general pattern, 
the calculation of Dn}-{c) amounts to sum the multiple series 

- hrrin 


D^nli.c) = —j- V] (mi ■ ■ ... ,m„) 

’ n\ 


mi! ■ ■ ■ m„! 


But we can also express it as an integral for a continuous measure. Indeed, 

d 


bn(t) = e ) e 


(54) 


(55) 


SO that (— l)”'e ^bn{—t) is the inverse Mellin transform of s"'r(s). We can then write 


PC+2CO 




s^a ^r(s)(is (c > 0). 


We can choose c = 1, and setting s = 1 + in, we obtain 

/•OO 




2na 


{1 + + iv)dv 


( 66 ) 


(67) 


whence the integral representation of the hyperdeterminants associated to c'^{a) = 
(-l)’"e““6„(-a), 

.X^7i(n-l)/2 


D'P(c'(a)) = 


(27ra) 


A^\v)Yla-^^^T{l + ivj)di 
j=i 


(58) 


Since Dn \c) is homogeneous of degree n in the c, and isobaric of weight kn{n — 1) 
(w.r.t. the weight function tn(c„) = n), we have 

DW(c'(a)) = e—Df(c(-a)) (59) 

and comparing both expressions, we obtain the identity 

/. n 

n(n+l) 


A2'^(n)n«““'r(l+2n,)dn,-= (-l)T(27ra)’^ ^ A^^m) 


(-a)”"i+-+™" 

mi! ■ ■ -m^! 


(60) 


i=l mi,...,mn>0 

which amounts to the evaluation of the integral (|58D by the residue theorem. Since 

-iv)^ 


T(1 + iv) = exp < —i(ln a + 7 )n + C(^)' 


m>2 


m 


(61) 


it is tempting to make the choice a = e~^^ in order to cancel the linear term in the 
exponential, and to get a curious identity involving on the left the values of the Riemann 
zeta function at the integers, and on the right Euler’s constant: 

Pra{-ix) 


[ A2^(a:)exp J VC(m) 

U>2 


m 


dxi ■ ■ ■ dxr. 


n(n — l) V—^ 

= (-l)^(27re-^)- A 


2k/ 


m 


3-7\miH- \-mn 


md ■ ■ ■ mj. 


(62) 
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Here, Pm(x) = power-sums. 

Let us now indicate an application of the above representation. Instead of considering 
the shifted hyperdeterminants which are the moments of the measure dfi'{x) = 

x^dp{x), one can replace by an arbitrary monic polynomial Q{x) of degree r. 
For a good choice of Q{x), the Hankel hyperdeterminants of the moments of 
dp.'{x) = Q{x)dp{x) may bear a simple relation to the original ones. 

In the case at hand, numerical experiments quickly suggest that such a simple relation 
occurs only with the choice 

Q{x) = {x)r = x{x — 1) ■ • • (x — r -I- 1) (63) 


and that we have then 

This choice amounts to replace our original sequence c„ = bn{a) by 


cl = feMfal = 


^s(r,j)6„+j(a) 
j=0 

where the s{r,j) are the Stirling numbers of the hrst kind. From (^), we have 




r^C+ioO 


T(s)5^s(r,j)(-s)"+Ms 

j=0 


pc+ioo 

2Tli 

(- 1 ) 

2711 


C—lOO 

r pc-\-ioo 


(—s)"a ®F(s)(—s)(—s — 1) • • • (—s — r + l)ds 
(—s)"'a“^F(r -f s)ds . 


Hence, 


DP(c") = 


(- 1 )' 

2771 

(-ly 


f*c-\-ioo 


A^^(s) a ®-’F(r -|- Sj)ds 


j=i 


nc-\-r+ioQ 

' c+r—200 


2vrz J 


i=i 


which is equivalent to 


(64) 


(65) 


( 66 ) 


(67) 


4. Miscellaneous examples 
4.I. Hilbert hyperdeterminants 

Another classical example of a Hankel determinant which can be evaluated in closed 
form is the Hilbert determinant 
1 


iHj -I 


= Dl'Kc) 


( 68 ) 


7J = 1 








Hankel hyperdeterminants 
where 


Cn 


x^dx. 


Jo 

Thus, ([T0|) gives immediately 

In the simplest case n = 2, the generating series (|53|) gives 


{2k + l){2k + 2)' 

Setting Cn{a) = (n + a + 1)“^, we obtain in the same way 

^f(c(a)) = 45n(l + a,l;fc) 
nl 

For a = r, we obtain the hyperdeterminants Zi)i^r(c). 
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(69) 


(70) 

(71) 


(72) 


4-2. A class of Hankel-Wronskians 

In [^, 4.12.3, one hnds the Hankel determinants Dn\c) associated to the sequence 


Cn = ^/(^) fit) = 


1 i / (73) 

1 — e* 

whose particular case x = 1 gives back one of the determinants computed by Lawden 
23|. To investigate the hyperdeterminantal analogues, we will hnd it convenient to 


make the substitution t ivr — t, and to assume at hrst that —x = iV is a positive 
integer. Up to a trivial sign, we can now take fit) = (1 + e*)^, and our sequence is 


'^ / -1 
Cn = — 1 + e 


t\N 


N 


k=0 


N' 


^kt 


N 


= (i+e‘rj2'‘' 


k=0 

N 


N 




k=0 


1 + e* 
p^(l-p)^-^ 


1 - 


1 + 


N-k 


(74) 


where p = e*/(1 + e*). That is, (c„) is the moment sequence of the binomial distribution, 
for which the Krawtchouk polynomials Knix]p,N) are orthogonal. We have 

N 


E 

k=0 


p>^il-p)^-^Kmik)K^ik) = 


(—l)"'n! — p 

i-N)n 


(75) 


whilst the monic polynomials Kn are related to the standard ones by 
/f„(i) = p”(-iV)„A'„(i). 


( 76 ) 
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Hence, 

\\k^f = {-lYn\p^{l-pY{-N\ 
which gives for the Hankel determinant 

n—1 

= (1 


13 


( 77 ) 


3=0 

t\n{n-l)/2 ^-1 


-e*) 


+ gt^n(-W+n-l) J.J 




(78) 


which agrees with Theorem 4.59 of 1^^ after substituting back t ^ iir — t and N = —x, 
and can be extended as usual to values of x ranging over the whole complex plane by 
means of Carlson’s theorem. 

The intepretation in terms of the Krawtchouk polynomials allows one to go one step 
further and to hnd a closed form for the D^^\. Indeed, we know that 

D'il = det(X„)D(') (79) 

where is the operator of multiplication by x followed by the orthogonal projection 
on the subspace spanned by the hrst n Krawtchouk polynomials. The matrix of Xn can 
be read directly on the three-term recurrence 

xkn{x) = kn+i{x) + \p{N -n)+ n(l - p)]kr,{x) + np{l -p){N + 1- n)kn-i{x) (80) 
which yields the tridiagonal matrix 

p 1 0 0 ■■■ 0 0 ^ 

a + p 1 0 ■■■ 0 0 

0 2A(/i-2) 2a + P 1 ■■■ 0 0 




(81) 


0 0 0 {n — l)X{ja — n + 1) {n — l)a + P _ 

where a = 1 — 2p, P = Np, X = p{l — p) and pi = N + 1. The three-term recurrence 
for the tridiagonal determinants is easily solved by means of a generating function, and 
one hnds 

det(X0 = (82) 

The other shifted determinants Dnl- can in principle be calculated by the same method, 
but is does not seem possible to solve the recurrences in closed form, and indeed, 
numerical calculations show that no nice factorised expression can be expected, except 
in the special case r = 2 and N = — 1, which gives back another one of Lawden’s 
determinants 1^. For the operator Xn \ multiplication by followed by projection. 


we obtain, for = — 1, 


det(X^2)) _ 




- 1 ) 


n+1 


(1 + 

but no other case seems to lead to an interesting formula. But this does not rule out 
the possibility of replacing x^ by another monic polynomial Q{x) of degree r. 
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To investigate this possibility, we shall adopt the same strategy as in the case of Bell 
polynomials, and look for an integral representation of onr sequence. Once again, it will 
be convenient to work with a slightly modihed (but equivalent) sequence 

an = ^9{t) where g{t) = (1 - e"*)^ = /(ivr - t). (84) 

The Laplace transform of g{t) is 


G{s)= e-^\l-e-Tdt= u^-\l - u)^du 


= B(s,iV + l) = 


N\ 


Hence, 


and 


rc-\-ioo 


s{s + 1) ■ ■ ■ {s + N) ' 


^n. — 


27ri 


s(s + 1) ■ ■ • (s + N) 


ds 




r‘C +200 


n\{2Triy 


A“(s) n 




=7 ^ji^j + 1 ) ■ ■ ■ {sj + N) 


Remark that by ([7^), we know that Dn-^r{a) is also equal to the hnite sum 

kr. 


^ \ nj 

which is indeed the value of the integral (|87|) according to the residue theorem. 


From (|87|) , we see that if we replace (a^) by the new sequence 

r-c+ioo 


1 pc-\-ioo 

^ / s-Q{s) 

ZTTZ J c—ioQ 


s(s + 1) ■ ■ ■ (s + N) 


ds 


where Q{s) = the hyperdeterminants 
d<P(r,N) = DP(a') 


(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 


satisfy to 

/ ]\T\ ^ 

40 (r; AT) = ( j N - r) (91) 

that is, are expressible in terms of the unshifted hyperdeterminants of the original 
sequence, with parameter N — r. 

As in the case of Bell polynomials, the Dn\c) can be reduced to a Pfaffian. 


5. Examples involving orthogonal polynomials 


Hankel determinants associated to sequences of the form c„ = Qn(x), where (Qn) is a 
family of orthogonal polynomials, have been called Turanians by Karlin and Szego, who 
computed their values for the classical families |^. Recent references on this subject 
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can be found in [Q, where these results have been generalized by a different method 
based on a little-known determinantal identity due to Turnbull. 

In this section, we will calculate the hyperdeterminantal analogues of the Turanians 
evaluated in 0 - As in the preceding section, we will make use of the integral 
representations of the classical orthogonal polynomials. Interestingly enough, Selberg’s 
formula will not be sufficient to deal with these cases, and we will have to rely upon one 


of its extensions, which is due to Kaneko |I5 


5.1. Kaneko’s integral and its variants 

The required integral formula involves the generalized Jacobi polynomials p^'^’'^{y) 
IjT^, ^ |, | 2 ^, which are the symmetric polynomials in r variables (yi ,... ,yr) obtained 
by applying the Gram-Schmidt process to the basis of monomial symmetric functions 
m^{y) (ordered by the condition p < n ii\p\ < |zz|,or \p\ = |z/| and p precedes z/ for the 
reverse lexicographic order) with respect to the measure 

r 

dp^'^^iy) = \K{y)\^^^^W(l-yiY{l + yifdyi---dyr (92) 

i=l 

on [—1,1]’’, normalized by the condition that the leading term of is m„(y). Let 

(93) 


i=i j=i 

Kaneko’s formula reads 


« n 

hoy- i=i 


\K{x)y‘^dxi ■ ■ ■ dxr, 


= 2 " 


S„(a, 6 ; c)p“;y (1 - 2 j/,,..., 1 - 2 ^,) 


(94) 


where a = ^ — 1 , /3 = ^ — 1 and 7 = c — 

The multivariate Jacobi polynomials indexed by rectangular partitions can be expressed 
in terms of generalized hypergeometric functions. This expression is simpler in terms of 
the polynomials 

(95) 

which are orthogonal on [0,1]'’ for the Selberg measure with parameters (a -|- 1, 6 -|- 1, c). 
For a rectangular partition k = (n’’), 

71 /c) f -n; a + b + s + n 
a + s 


-7/ • I'l _ Tpjy) ( 

^{n^) 1^1’ • • • 1 2/r, ^) — 2^1 I 


2/1 


y ■ ■ ■ iVr 


(96) 


where s = 1 -|- (r — l)c. The generalized hypergeometric functions associated with Jack 
polynomials CiJ\yi,... ,yr) are defined by 0 0 


p(a) ( 0-1 ■ ■ - Clp 


1 [ai]i ^ 

yi,...,yr] =2^-i2^ 


i(«) 


^ n! ^ ■ ■ - Ih F' 

n>0 \K\=n 


l(“) 




(97) 
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where 


1{k) 


HL”' = 


2=1 




We note that in the case n = ijf), the denominator of ( pS]) is given by (|9lD as 

,S'n(a ,6 + r, c) 


a,o,c— i 


Pk 




16 


(98) 


(99) 


Sn{a,b,c) 

This formula is needed to calculate the degenerate cases of Kaneko’s integral 
corresponding to Laguerre and Hermite polynomials. 

The generalized Laguerre polynomials (?/;«) are dehned by (see also |^) 


K(vu a) = Jim (f. • • •. f i ») 


( 100 ) 


(we use there the convention of [^]). Let LSn{a, c) denote the Laguerre-Selberg integral 
(p5|). One can deduce from (pi]) the Laguerre version of Kaneko’s integral. Indeed, 
Kaneko’s formula can also be written as 0 


'[ 0 , 1 ]’ 


R{x,y)\A{x)\'^''Yl^i ^dxi 


2=1 


= Sn{a + r, 6, c)2Kf 


(c) 


n; -(a -|- 6 -|-r — l)-|-n — 1 




lia + r) - 1 

Setting Xi = Ui/L, yi = Vi/L and letting L —> cxo in this formula we obtain 


( 101 ) 


' (a,cx))’ 


i?(M,n)|A(M)p^ ^"dui 


2 = 1 


= lim inr+cn{n-l)+n(a-l)+ng^^^ + r, L + 1, c) X 


L—*oo 


(c) f -n,l{a + r + L)+n-l 


From (1^) and (|10(]|) we have, setting b' = — 1, 


L''"' L 


( 102 ) 


lim 2^1 

L—^■oo 


F) ( + r + L) + n - I 


i(a-|-r- 1 ) 


Vi Vr 

T''"'T 

/ I j,/ I „/ 


= lim 2 F\ 

6 '—>O0 


= LiC-r--7-7c) 

c c 


(c) f — 1 , -|- n 


a' + s' 


Vl 


cb' + c — 1 cb' -\- c — 1 


where a' = ^ - 1 and s' = 1 + ^. 

C C 

On another hand, 

lim inr+cn{n-l)+n{a-l)+ng^^^ + ^ + 1 - c) = LSn{a + r, c) 


L—>oo 


(103) 


(104) 


Finally, we get 


'( 0 ,oo)’ 


R{x, y)A^^{x) '"'dxi = LS'„(a r, 


fyi yr 


2 = 1 


^ ^ c c 


(105) 
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Similarly, an appropriate limit of (Q) yields 
[ R{x, y)A^^{x) e-^Uxi = ^ 


i=l 


i=i 


X 


where the generalized Hermite polynomials H^iy^a) are dehned by 
H^{yi, • ■ ■, a) = lim (-A/^)l^lL“(a + yiV^, ■ • •, a + yr\^] a). 


( 106 ) 

( 107 ) 


We follow here the convention of [Q]. 

Kaneko’s identity can be interpreted as a generalization of Heine’s integral 
representation of orthogonal polynomials in the Jacobi case. Indeed, it can be rewritten 
as 

/ \A{x)\‘^^d^it{xi)---dnt{xn) ( 108 ) 

where d/atix) = - a;)(l - 3 :)“-^! + xy-\ = 2cn(n-i)+n(a+6+r-2)^^('^^ 

a = c{a + 1), b = c{P + 1 ), c = 7 + 1 . 

Hence, when c = fc is a positive integer, the symmetric Jacobi polynomials indexed by 
rectangular partitions are expressible as hyperdeterminants 



= n\D>P(c(t)) 


( 109 ) 


where 

Cm{t) = j x^dutix). ( 110 ) 

This can be regarded as a generalization of the classical determinantal expression of the 
orthogonal polynomials in terms of the moments. 

The extension of these identities to non-rectangular partitions or to other measures 
appears to be unknown. 


5.2. The case k = 1: Leclerc’s identity 

On another hand, for k = 1, Kaneko’s representation can be extended to general 
orthogonal polynomials. Let /i be any linear functional such that the bilinear form 
if yd) = Tifd) is non degenerate, and extend it as above to functions of n variables 
xi,...,Xn by setting p.n{fi{xi) ■ ■ ■ fn{xn)) = /i(/i) ■ Let pf\xi, ...,Xn) be the 

basis of symmetric polynomials obtained by applying the Gram-Schmidt process to the 
monomial basis with respect to the scalar product 

(/, 9)k = hn {A‘^’'{x)f{x)g{x)) (111) 

with leading term m^. In this section, we shall use the representation of partitions by 
weakly increasing sequences A = (0 < Ai < ... < Xn) instead of the usual one (this will 
be more convenient for the indexing of minors). The one-variable polynomials Pm(x) 
are the monic orthogonal polynomials associated with p. 
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When k = 1, one has 

( 1 ) / \ 

Px [x) = 


Dx{x) 

A{x) 


( 112 ) 


where the alternants (Slater determinants) D\{x) = det(pAi+i-i(a^j)) form the natural 
basis of antisymmetric orthogonal polynomials for Now, we can write 

- n n(fe - ^i) - 

The analogue of Kaneko’s integral in this context is the scalar product 

x)) = (1, R{y, x))i. (114) 

Expressing A{x,y) in terms of the one-variable monic orthogonal polynomials Pm as 

A{x,y) = det{pi_i{xj)\pi_i{yk)) (115) 

and taking the Laplace expansion of this determinant of order n + r with respect to its 
hrst n columns (containing the variables Xj), we hnd that 

= At ^ = 5^(-l)'“W(a^)P/3(2/) 

where the sum runs over all pairs of partitions 

a = (0 < ai < ... < Q!„), 13 = {d < j3i< ... < j3r) 


(116) 


(117) 


such that (tti + 1,..., an + n, Pi + ^ Pr + r) is a permutation of (1, 2 ,..., n -|- r), 

in which case (—l)l“l is its sign. In particular, for a = 0, is the rectangular partition 
P = {rf)-, so that 


{l,R{y,x))i = 5^(-l)'“'pJ^^(2/)tf 

a,P 


(118) 


This equation contains as a special case Theorem 1 of |^, which in turn implies all 
the identities of Karlin and Szego as well as many other ones. Indeed, taking the limit 
?/j —*• M, i = 1,..., r in (|112|) , we obtain a Wronskian of one-variable polynomials 

W{PXi,PX2+U ■ ■ ■ ,PXr+r-l){u) 


( 1 ) / \ 


(cf. 7.1.1 p. 107), so that 

Hn ( A^(x) ^(y - XiY j = /in(A^) 


l!2 !---(r- 1)! 

W{pn,...,Pn+r-l){y) 


i=l 


l!2!---(r-l)! 


But we have also 

Hn I A‘^{x) TT(|/ - XiY 1 = (-1)”''^! det(c,.+i+j(?/))| 


n—1 

0 


(119) 


( 120 ) 


( 121 ) 
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where 


Cm( 2 /) = vT] = ^ ^. 

3=0 


m 


( 122 ) 


and (|121|) is precisely Theorem 1 of p6 


The equality of the right-hand sides of 
Applying this identity (with ?/ = 0) to the case where the moments c„ are the Bell 
polynomials hn{a) = so that Pn{x) = Cn\x), we obtain (^). 

This suggest the conjecture that in general pn{^‘^’‘{x)R{y,x)) should be expressible as 
where the ^ are the symmetric orthogonal polynomials for another 
functional fi' related to qt in some natural way. 

5. 3. Hyperturdnians of Legendre polynomials 

Let us start, as in 0 , with the Legendre polynomials Pn{x). Laplace’s integral 
representation 

Pn{x) = — [ {x + coscfVx'^ — l)"'d(j) (123) 

^ Jo 

can be rewritten as 

Pn{x)= [ Pdp{t) (124) 

J a 

where a = x — \/x‘^ — 1, h = x + ^x^ — 1, and dplt) = 7 r“^(f — a )“^/^(6 — Hence, 

1 S’ 

Df\P(x))^— (125) 

which under the substitution ti = {b — a)ui -|- a reduces to the Selberg integral 

1 n n 

D^n'(no) = -r^(i> - / A“(«) n“r‘'''(l - n0-''^du, 


niTT'' 




Now, the shifted polynomials c„ = Pr+n{x) are the moments 

Pr+n{x) = f Pdprit) 

J a 

where diiy{t) — TT~^V'[t — a)~^P(b — t)~^Pdt, so that 


(126) 


(127) 


DW(P(a)) = 


Jla,b] 

_ ^'jkn{n—l)+rn 


nlTT^' 


A^^(f)(ti.. .4)'']^(fi - a) ^^^(b-ti) ^I'^dti 

i=\ 

n 

A^^(m) ]^(Mi + vYu~^^‘^{l - Ui)~^^‘^dui 


(128) 


'[ 0 , 1 ]" 


2=1 
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where v = This is of the form (^) with yi = y 2 = ■ ■ ■ = Vr = —v, whence, since 


mvr" 


X P, 


«,/3,7 


X 


X 


(»') 


(129) 


'X^ — 1 — 1, 

where Q! = ^ — 1 . = and 'j = k — 

It is instructive to have a look at the case k = 1. Here, a = (5 = —\ and the generalized 
Jacobi polynomials are the symmetric orthogonal polynomials for the measure 

dyi 


dfi{y) = A'^iy) JJ 


(130) 


^=l Vi - yf 

For r = 1, the orthogonal polynomials are the Chebyshev polynomials Tn{y), and 
the symmetric orthogonal polynomials for (|130|) are the D^{y)/A{y) formed from 
corresponding monic polynomials t^- Taking the limit of these expressions for 
{yi,...,yr) —>■ ('Ci-'-iOi where ^ we obtain a Wronskian of Chebyshev 

polynomials evaluated at which is precisely the expression of the Turanian found 
by Karlin and Szego (see also p 6 [|). 


5.4- Laguerre 

We start with the hypergeometric representation of the monic Laguerre polynomials 

-n, b 




-n 


Qj 


X ) = lim 2^1 . 

6—^00 V a + 1 


(131) 


The second part of this equality leads to write each shifted hyperturanian as the limit 
of a Kaneko integral, which gives after simplihcation 


d£(lI“>) = 


lim 


kn{n—l)-^nr n—1 


n\k\^ 6-100 \ b J 

A_1 a-fc+l I b 

' ' ^ h ' 


n 

i=o 


(6)jfc+r(a - & + l)jfc(jfc + k)\ 


(ci T 


X 


X Pd 


X X 


(132) 


From (1^) , we see that this can be written as a generalized hypergeometric function 


r)(k)(T{a)\ _ ( +nr knin-l)+nr ^ TT 


{jk + k)\ 


X 




(133) 


In particular, if r = 0, we obtain 


D^(L(“)) = 1 ) ^ 


n-l 


_TT ijk + k)\ 

n\k\^ (a + l)fc(n+i-i) 


(134) 
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We start as above with the representation of the monic Hermite polynomials as limits 
of hypergeometric functions 


HJx) = lim 0 2 2 F 1 


—n, 2a 


lfi-4 

2 V ^/i 


(135) 


We can then write the shifted hyperturanian as the limit of a Kaneko integral. 
If r > 0, from (|106|), one hnds 


1 _ 


i=i 


In the simplest case r = 0, we obtain 


~ 1 1 ^ 


(137) 


J = 1 


(\ iV 

{N = which can be found in (17.6.9). Indeed, expanding A{xY^, one has 


. , dx'^, 

Vj = l ^ . 


N 




a'l,---,o'2fce6n 


X 


X 




N 


d'^h 


11:^ If 


X, 


CTi (i)H- \-a-2k{i)-2k 


Jl - ■ - In J dx^j 

But for each monomial x\^ ■ ■ ■ x^ appearing in the previous formula, we obtain 

N 


(138) 


E 


h---lnj ,J[dx 


nliin 


x: 


.pi 


Xi=0 


N\ nr=i In if Pi i® 

2 • 

0 otherwise. 


It follows that 




N 


^ = i-lf 7112 “^^N\dI^\H) 


. , dx‘^, 

,j=i j 


Hd^)! 


k\’ 


(139) 

(140) 

(141) 


J = 1 


5.6. Charlier 

The monic Charlier polynomials C^\x) = n\Ll^^~^\a) are given by the exponential 
generating function 

j-n 

Ect‘’(^);^ = «'“‘(i + i)'- (1^2) 

n>0 
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(143) 


r(-x) Jo 

From this, we get easily the Hankel hyperdeterminants associated to the sequence 
Cn = Ci°'\x). The result can be cast in the form 


n—1 




(k + kj)! r(—x + kj) 


j=0 


k!r(-x) 


In particular, for the determinants. 


n—1 


c<‘>(c)= n 


j-o 


r(j + i)r(j -x) 

r(-^) 


(144) 


(145) 


The shifted hyperdetermiiiants can be similarly evaluated via Kaneko’s identity 

in the Laguerre form, 

n—1 




(146) 


j=0 


For the shifted determinants, we get a Wronskian of Laguerre polynomials which is 
easily seen to be equivalent to the evaluation given in [|^ . 


5.7. Meixner 

The Meixner polynomials admit the integral representation (|^ p. 448) 


0„(-x;/3,7) = 


r(/3) 


t^-\l-tf 


—x—1 


1 + ( --1 

7 


dt. 


T{f3-x)T{x) Jq 

From this, one deduces the hyperdeterminants associated to Cn = (3,^) 

n(fc) ( ] = L ( kzi F(x + jk)T{f3 -x + jk)T{(3)T{jk + k + 1) 

n! V 7 y f}^T{x)T{p-x)T{P + {n + j-l)k)T{k + l)' 

Kaneko’s identity gives directly the shifted hyperdeterminants as 


(147) 


(148) 


nW( 1 = 1 /IpTr'"”’"*”’’ ff T(x+jk + r)T{l3- x+jk)T(p)T(jk + k+l) 
n! f 7 ) f|r(j-)r(/3-a)r((3 + r+(n + j-l)k)r(k+l) 




X 


^ P(nn(Pr--,P', k) 


(149) 


with a = I — 1, 6 = — 1 and p = 

5.8. Krawtchouk 

The Krawtchouk polynomials are given by 

Kn{x-,p,N) = 2 -Fi(-n, -X] -N; -). 

P 


( 150 ) 
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Assuming at first that —N is not a negative integer, we can write down an integral 
representation 


Kn{x;p,N) = 


T{-N) 


f 


T{-x)T{-N + x) Jo 
which leads immediately to the evaluation 

1 / ^i-N) 


-x-in 4^-v+x-l I i_l] dt 


(i-ty 


dP(k) = 


n! \T{—x)T{—N + x) 

After simplihcation, we hnd the expression 


^ \ kn(n—l) 

-) Sn{-x,x - N,k). 

pj 


dP(K) = 


^]pkn{n—l) 


n—1 

n 

j=0 


'-x)jk{N - x)jk{jk + k)\ 

(A^)A:(n+j —l)+r^' 


(151) 


(162) 


(163) 


which is well dehned for integral provided that all the elements of the 
hyperdeterminant are also dehned (recall that Kn is dehned only for n = 0 ,..., A^). 
The shifted hyperturanians can be evaluated from Kaneko’s integral, 

(-ir 


Yir /i \ kn{n—l)-\-nr 

^ C-) 

X (-a^) jfc+r(-A' + x)jk{kj + k)\ p{a,l3) 
j=0 


pp^>{p,...,p-,k) 


(154) 


( A^)fc(n+j —l)+r 

■ 1. In particular, for fc = 1, we obtain a Wronskian 


where a = —| — 1 and (3 = ^ 

of Jacobi polynomials with parameters —x — 1 and —N + x — 1. 


6. Hankel hyper determinants and symmetric functions 


In this section, we shall give an expression of the Hankel hyperdeterminant Dn\c) in 
terms of symmetric functions. Precisely, we suppose here that = hn(x), the n-th 
complete homogeneous symmetric function of some auxiliary set of variables x = {x*}, 
and our aim is to obtain an expression of the symmetric function Dn \h) in terms of the 
Schur functions sa(x) (see for notation). It turns out that this problem is equivalent 
to hnding the Schur expansion of the even powers of the Vandermonde determinant, a 
difficult problem which has been thoroughly discussed in recent years, mainly in view 
of its potential applications to Laughlin’s theory of the fractional quantum Hall effect 
(see and references therein). 


Since Dn\h) is a homogeneous polynomial of degree n in the hi, its Schur expansion 
will only involve partitions of length at most n. We can therefore assume that 

X = {Xi, • • • ,Xn}. 

It will be convenient to work with Laurent polynomials in x. In particular, for each 
vector A = (Ai, • ■ ■, A„) G Z"", we dehne the augmented monomial symmetric function 

„crA 


m\ = 2_^ X 
cre 6 „ 


where (t(Ai, • • •, A„) = (A^(i), • • •, A^(„)) and x^ = x^ 
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Let (f) be the linear map sending rhx to if A G N"’ and 0 otherwise. As the set of 
the rhx, with A a decreasing sequence, is a basis of the space of symmetric Laurent 
polynomials the map 0 is well dehned. 

Let us consider now the alternants 


ax = 




(TG&n 


The image by (j) of the product of 2k alternants is a hyperdeterminant. Since 

2k 


i =0 


O'!, 


,(T2k-l&&n 


we get 



= Det 


2k 


I \( 2 *) 


(165) 


(156) 


(157) 

It 


The case where fc = 1 is well know and can be found as an exercise in the book 
is shown there that for any symmetric function /, we have 

(j){fasa-s) = f , (158) 

where 5 = (n — 1,..., 2,1, 0). In particular, 

Dlf\h) = Det2t(5„+...+i,.)r‘ = M) = ■('(A(a^)"‘) 


2k 

nin — l) " 1 —r 

[-1)^ l[xr ^Aix] 

i=l 


2{k-l) 


astt-s 


= {-i)^l[xrA(xr-'K 


(169) 


2=1 


Since we are working with n variables, the effect of the factor nr=i f® shift the 

parts of the partitions occuring in the Schur expansion of by n — 1. 

The Schur expansions of A^, which determines all Hankel hyperdeterminants of order 
4, have been computed up to 9 variables in and recently up to 10 variables by 
Wybourne. Using the Littlewood-Richardson rule, it is then possible to compute the 
powers A^^ for small values of k. The hrst cases are 

^2 \^) — “Ssi + 3S2 


’22 


('OA 

-^3 W — ~'S642 + 3S633 + 3S552 — 6 S 543 + I 5 S 444 

Df\h) = -S 51 + 5s 42 - 10s, 


33 


D 


(3)i 


3 (h) — —Si062 + 5 Sio53 ~ rusi044 “T OS 972 


10 sin 44 + 5S972 — 20S963 + 25S954 
10S882 + 25S873 + 15S864 ~ lOOSsss — IOOS 774 + I 6 OS 765 


280s, 


666 


D 2 — —Syi T 7sg2 — 2 IS 53 + 35s 


44 


In terms of the elementary symmetric functions e„ and the power sums pn, this identity 
can be rewritten as 


DP(h) = (-l)"'”-‘''^er'‘det(p„_.+,) 


fe-1 


(160) 
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As an illustration of ( |159| ), let us consider the case where c„ = Un{x), the nth Chebyshev 
polynomial of the second kind. From the generating function 

1 


J2Un{x)t^ 


n>0 


1 — 2xt + 


we see that Un{x) = hn{xi,X 2 ), where Xi + X 2 = 2x and XiX 2 = 1. Hence, 
= {xi + X 2 )‘^ — 4x1X2 = 4(x^ — 1), so that = — [4(x^ — Comparing 

with (0), we obtain the identity 




i=o 


0 


otherwise. 


( 161 ) 


Specializing to the Fibonacci numbers /„ = {—i)^Un{i/2), (normalized such that 
fo = fi = 1), we hnd 

, 2k 


j=0 


2k 




'fc-i 


( 162 ) 


6.1. An application: inverse factorials 

The symmetric function approach allows us to handle the case 

Cn = —. (163) 

n! 

Indeed, the image of a Schur function s\ under the specialisation h->■ c„ is equal to 

the scalar product 
1 


N\ 




(164) 


where A = |A|. If A has at most n parts, we can interpret the above as the scalar 
product of rational GL(n)-characters, dehned by 


{f{x),g{x))GL{n) = 


n 


■CT|/(x)n(l-fi)9(*)| 


(165) 


where CT means the constant term, and x = (x^^,..., x“^). 
Hence, we can write 

1 


l(k,n) = Bf(c) = _ ((-1)=^(^,. .. a:„)"-'A“-2(a:), (a:, + ■ ■ • + 


[kn{n — 1)] 

. kn(n — l) , ^ 

(cT.„)‘<"-')n(i-^) 


[kn{n — 1)]! n! 




Xi 


Xr. 


kn{n—l) 


.. kn(n — l) 

- 1 ) ■ 


((ll . . . x„f , (xi + ... + 


[kn{n — 1)]! 

where in the last equation we have introduced Macdonald’s second scalar product { , )'„ 
associated to Jack polynomials in n variables, with parameter a = ^ (see |2^, (10.35) 
p. 183). 


(166) 
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Now, (Xi • = P(k(n- 1 ))^^ 


{k{n-l)y 

(xi + ■ ■ ■ + ci^\ 


X) 


K.\-kn{n—l) 


where, if k h A^, 


CkW 


(167) 


(168) 


with 


C/7(tt) = JJ (tt(«7 - j) + K'j-i + l). (169) 

{i,j)£K 

Therefore, if we denote by u the rectangular partition {k{n — 1))” of weight N = 
kn{n — 1), the Hankel hyperdeterminant is given by 

- 1 )^ 


(170) 


which can be evaluated in closed form thanks to equation (10.37) of p. 183. This 
yields 

[-lfP^-^^P{kn)\ {ki)\ 

{k{n + i-l))\ 


l{k,n) = 


n\{k\y 


(171) 


The case k = 1 could have been obtained in a simpler way from the hook-length formula 
giving the dimensions of the irreducible representations of the symmetric group. 


7. Conclusion 

We have demonstrated that the calculation of Hankel hyperdeterminants amounts to the 
evaluation of an interesting class of multidimensional integrals, including Selberg’s and 
Kaneko’s ones, and more generally, the partition functions of one-dimensional Coulomb 
systems with logarithmic potential. We have presented a series of examples which can 
be evaluated more or less directly from known results, and obtained on our way a 
unihed presentation of many Hankel determinants, including some new cases. However, 
obtaining new integrals from algebraic or combinatorial evaluation of hyperdeterminants 
would be more interesting. A few examples are presented in this paper, and we expect 
more from a systematic study of Hankel hyperdeterminants from an invariant theory 
point of view. Indeed, any generalization of one of the various tricks working with 
ordinary Hankel determinants would immediately lead to new interesting integrals. 


Appendix A. Hyperdeterminantal aspects of Selberg’s original proof 

The Selberg integral can be deduced from the Hilbert, factorial or inverse factorial 
hyperdeterminants. Actually, the evaluation of any non-trivial class of Hankel 
hyperdeterminants would lead either to Selberg’s intergral in full generality, or to some 
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interesting generalization. This is already apparent in Selberg’s original proof. The first 
part of his proof can be translated in terms of hyperdeterminants in the following way. 
First, we write Selberg’s integral, for c = k an integer, as a hyperdeterminant 

S'„(a, b, k) = n\B{a, {3)"'D^\c{a, a + j3)) (A.l) 

where c„(a, h) = and {a)n = a{a + 1) ■ ■ ■ (a + n — 1) is the Pochhammer symbol. 

Computing S'„(a, b, k) amounts to find a closed form for Dn\c{a, b)). 

Expanding the hyperdeterminant, we find an expression of the type 


Dl^\c{a,b)) = 




n 


(®)fc(m—1) 




(®)im ib)k{n+m-2) 


(A.2) 


(A.3) 


(P)k{n+m-2) j i(^)k{m-l){b)j^ 

where J runs over the integer vectors J = such that ji < j 2 < ''' jni 

ii + ■ ■ ■ + in = kn{n — 1), and cj G Z. 

A straightforward investigation of the 2A;-uplets of permutations giving a non zero cj 
implies that for each m G {1, ...,n}, one has 

k{m — 1) < im < k{n + m — 2). (A.4) 

Hence, our hyperdeterminant can be written as a product 


DfWa,4))= n 


{,oi)k(m—i)ib (k)k(m—i) P(^a,b a) 

— X 


m=l 


k(n-\-m—2) 


Qib - a) 


(A.5) 


where P(a, b) is a polynomial whose degree in b is at most and Q{b) is a 

polynomial of degree 

Since Dn \c{a, b)) is symmetric in a and & — a, we see that the ratio = Q,{n, k) 

is independent of a and b, so that 


Bf(c(a,6)) = aKl')fl^%^'’ “I"'”-* 

™=i 


(n+m—2) 


(A.6) 


(A.7) 


Now, setting a = 1 and b = 2 we obtain 

{k{m-l)y/{b)kim+n-2) 

that is, Selberg’s integral can be deduced in full generality from the Hilbert 
hyperdeterminant. 

Also, observing that the case of inverse factorials is given by the limit 
l{k,n) = lim L-’^^^^-^^Diy\c{L + l,l)) 


L—>oo 


we have 


fl<‘>(c(a.6))=(-i)^nn 


{(^)k{ m— i){b-a)ki m—1) 

•jl (Hn + m- 2))\{b)k(n+m-2) 


l{k, n) 


(A,8) 


(A.9) 
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In the same way, the Hankel hyperdeterminant of factorial numbers 
¥{n,k) = lim + 1)) 

L—^■oo 

gives another equivalent identity 
Df(c(a,h)) = n 


(A.IO) 


(A.ll) 


Appendix B. Possible generalizations 

More generally, if we start with a hypergeometric moment sequence c„ = where 

P{n) = and Q{n) = polynomials in n, a similar analysis 

leads to an expression of the hyperdeterminant Dn \ c ) as a product 


km —1 


n —1 


n 


m=0 


(B.l) 


n 


j =0 


where Rn'^ (a; 6) is a polynomial of degree at most in both sets of variables 

a = {oo, • ■ ■, ttr} and b = {Bq, ■ ■ ■ bg} and whose coefficients are in Z. 

In the most general case, R^\a',b) cannot be factorized, and seems difficult to compute. 
However, in some simple cases, we can give a closed form. 

Suppose that P{n) = ai? + bn + a and (^(n) = 1, we hnd 

2fc+l 

(B.2) 

i = k -\-2 


Rf\a,b,c]l) = JJ {b + cj). 


Let us give now two examples involving combinatorial numbers. 

(3) ( ) 

The tri-Catalan numbers Cn admit a representation as moments 

gSn+l 2^ 


32 


cp) = 


27 

4 


2^/371 
2 , 


2n -I- 1 


x^djjL{x) 


(B.3) 


where duix) = v^ 2 f 2?(27+3V8i 123 :)? dx. It follows that our hvDerdeterminant has 

’ 1271 x3(27+3V81-12x)3 

the integral representation 




4 


'0 




(B.4) 


i=l 


which looks rather difficult to compute. Nevertheless, our previous remarks allows to 
start the calculation 

n T-rfcm-l/_2 , 77 ■ 10 -2 , g.' 3 ^ 

Bf(C®) = n ^ 1 ■R!.'‘’(-2.11.-18.9;0.-1.4) (B.6) 

m=0 ffi=0 d) 

and it remains to hnd a closed form for i?P^(ai, 02 , 03 , 04 ; &i, & 2 , ^s)- When k = 1, the 


result is known [3S . 
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If Cn = {2n)\ we have 




JO v-^* 

(B.6) 


n km—1 


_ 2^kn{n- 

-b Himk-iy H {2j - l)Riy\0, -2,4; 0) 

(B.7) 


m=0 j=0 


Let us give now some polynomials Iin\a,b) for various values of n, fc, a and b: 
R^\ao, oi, 02; 60, bi) = —aobi + aibo + 80260 + 20261 

R ^\ ao , oi, 02; 60,61) = 6(—2ooOi6o6i — 1009026061 + 1501026061 — 090161 — 1509026^ 
-|-9(2ifZ2&Q H“ ^^0^1 H“ ^1^0 20q-2^q “h 24 o- 2&^ “h “1“ 50fZ2^o^l) 

^1: ^2: ^3? ) ~ “h 9 QiQ- 2 “t“ 200*2 ^0^3 “I” 350i03 “h 150o2^3 ~t“ 2740^) 

i?3^\ao, oi, 02;) = 16(94251 ooOi 02 + 5525ooojo2 + 48 oqOi 02 + 1853066oooi02 

-I-6O3IOI00O2O2 25518oqOi02 7123ooo^02 522oqO|02 32783900^02 

+15303958o?o^ + 384ogoi + 41544o^o^ + 211ooo^ + lOogo^ + 592o? 
+371151360^ + 21786960^00 + 37277876o^oi + 385834ofo^ + 23654o^02) 

(; bo , 6i, 62,63) = -6(6063 - 6? - II6162 - 456163 - 3061 - 2506263 - 52461) 


Appendix C. Pseudo-hyperdeterminants 


For tensors of odd order, another notion of hyperdeterminant is considered for example 
in [03 


Det^ ()o ^ ^ ■ ■ ■ *^2fc) n ^icri(i)---(T2fe(*)' (C.l) 

(Tl---(T2fc i = l 

Note that this polynomial has not the same invariance properties as the 
hyperdeterminant under linear transformations. 

When = Cmij+i 2 +---+i 2 fe+i) being the moments of a mesure dfi{x), this 

hyperdeterminant can be expressed a multiple integral involving an even power of the 
Vandermonde determinant 


. n—1 
'0 

n 


mi +( Ti { i )+--- o - 2k { i )- 2k 


-^rni^) Det_|_ (Cmij^+i 2 H-|-* 2 fc) 

It n 

= n / 

CTl--t72fe *=1 “ 

pb 

= /•••/ 4^1 ■ ■ ■ (^ 2k ) n 

tJ CL J CL — — : _1 

( 7 \"- G2k ^ —-L 

. n 

d^,{xi) ■ ■ ■ dfi{xn) 


f*b rb ^ 


dfi{x) 


mi + ai { i ) + --- cr2k { i )- 2k 


djji^Xi) 


(C.2) 


'a J a 


2 = 1 


Obviously, one has 

■){k 
( 0 ") 


DW,(c)=n!Df(<;) 


(C.3) 
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and we can compnte other examples related to Selberg’s integral. The main tool is the 
system of differential eqnations verihed by the fnnctions = A(a;)^^ 11^=1 — 

XiY~^ (see Aomoto’s proof of a variant of the Selberg integral in |]^ or for example). 
Let ns list a few resnlts: 


(i) For Cn = one has 


T~\W 


»1 


Xf ■ ■ Xs^{xY^dxi ■ ■ ■ dx 


'0 ^0 


^ ^ 1 ^ 1 + {n-j)k (fc(l + j))!(fcj)!^ 

fc!” 2 + (2n — j — l)k (1 + (n + j — l)fc)! 

More generally, if c„ = one gets 

a + k{n — j) 


(ii) 




_^_n 

r(6 — a)^k\^ b + (2n — j — l)k 


X 


111 


n—1 

xn 

j=0 

If Cn = n!, one obtains 

^ m 


{k{l+j))\T{a + jk)T{b-a + jk) 
T{b + {n + j - l)k) 


n—1 


^(n,0"-)(c) = + Hn - j))Yl{H^ + j))KkjV- 


i=i 


i=o 


and 


+ -D(2m IS (c) J. k{2Tl TR j)) X 

m+s n—1 

i=i i=0 


(C.4) 


(C.5) 


(C.6) 


(C.7) 


References 

[1] K Aomoto, Jacobi Polynomials associated with Selberg Integrals, S.I.A.M. J. Math. Anal. 18 

(1987), 545-549. 

[2] T H Baker and P J Forrester, The Calogero-Sutherland model and generalized classical 

polynomials, Comm. Math. Phys. 188 (1997), no. 1, 175-216. 

[3] D Barsky, M Carpentier, Polynomes de Jacobi generalises et integrales de Selberg, The Electronic 

J. of Combinatorics 3 (2) (1996). 

[4] A I Barvinok, New algorithms for linear k-matroid intersection and matroid k parity problem. 

Math. Programming 69 (1995), 449-470. 

[5] A Cayley, On the theory of determinants, Trans. Cambridge Phil. Soc. VIII (1843), 1-16. 

[6] A Cayley, Memoire sur les hyperdeterminants, J. reine und angew. Math. 30 (1846), 1-37. 

[7] A Cayley, On the theory of per mutants, Cambridge and Dublin Math. J. VII (1851), 40-51. 

[8] A Debiard, Systeme differential hypergeometrigue et parties radiales des operateurs invariants 

des espaces symetrigues de type BCp, in Seminaire d’algebre, M-P Malliavin ed.. Springer 
Lect. Notes in Math. 1293 (1988), 42-124. 








Hankel hyperdeterminants 


31 


[9] M Desainte-Catherine and X G Viennot, Enumeration of certain Young tableaux with bounded 
heighf in Combinatoire enumerative (Montreal 1985), Springer Lect. Notes in Math. 1234 
(1986), 58-67. 

[10] P J Forrester, Random matrices, log-gases and the Calogero-Sutherland model, in Quantum many- 

body problems and representation theory, 97-181, MSJ Mem. 1, Math. Soc. Japan, Tokyo, 
1998. 

[11] I M Gelfand, M M Kapranov and A V Zelevinsky, Discriminants, Resultants and 

Multidimensional Determinants, Birkhauser, Boston, MA, 1994, 523 pp. 

[12] F Gherardelli, Osservazioni sugli iperdeterminanti, Istit. Lombardo Accad. Sci. Lett. Rend. A 

127 (1993), no. 1, 107-113 (1994) 

[13] P Haukkanen, Higher-dimensional GCD matrices. Linear Algebra and its Applications 170 

(1992), 53-63. 

[14] H E Heine, Handbuch der Kugelfunktionen, 2 vol., 2nd edition, Berlin, 1878, 1881. 

[15] J Kaneko, Selberg integrals and hypergeometric fuctions associated with Jack polynomials, S.I.A.M. 

J. Math. Analysis, 24 (1993), 1086-1110. 

[16] S Karlin, Total positivity, Stanford University Press, Stanford, 1968. 

[17] S Karlin and G Szego, On certain determinants whose elements are orthogonal polynomials, J. 

d’Analyse Math. 8 (1961), 1-157. 

[18] R Koekoek and R F Swarttouw, The Askey-scheme of hypergeometric polynomials and its q- 

analogue. Dept, of Technical Mathematics and Informatics, Report no. 98-17, Delft University 
of Technolgy, 1998. 

[19] A Koranyi, Hua-type integrals, hypergeometric functions and symmetric polynomials, in 

International symposium on number theory and analysis in memory of Hua Loo Keng, Gong, 
Sheng (ed.) et ah. Volume II: Analysis. Berlin: Springer- Verlag, 1991, 169-180. 

[20] C Krattenthaler, Advanced determinant calculus, in The Andrews Festschrift, Seminaire 

Lotharingien de Combinatoire 42 (1998), electronic, 67 pp. 

[21] J S Kung and G C Rota, The invariant theory of binary forms. Bull. Am. Math. Soc., New Ser. 

10 (1984), 27-85. 

[22] M Lassalle, Polyndmes de Jacobi generalises, C. R. Acad. Sci. Paris Ser. I 312 (1991), 425-428. 

[23] M Lassalle, Polyndmes de Laguerre generalises, C. R. Acad. Sci. Paris Ser. I 312 (1991), 725-728. 

[24] M Lassalle, Polyndmes de Hermite generalies, C. R. Acad. Sci. Paris Ser. I 313 (1991), 579-582. 

[25] D F Lawden, The function and associated polynomials, Proc. Cambridge Phil. Soc. 

47 (1951), 309-314. 

[26] B Leclerc, On certain formulas of Karlin and Szego, Seminaire Lotharingien de Combinatoire 41 

(1998), electronic, 21 pp. 

[27] J-G Luque and J-Y Thibon, Pfaffian and Hafnian identities in shuffle algebras, Adv. Appl. Math., 

to appear. 

[28] I G Macdonald, Symmetric functions and Hall polynomials, 2nd ed., Clarendon Press, Oxford, 

1995. 

[29] M L Metha, Matrix theory, Hindustan Pub. Corp. and Les Editions de Physique, 1989. 

[30] M L Mehta, Random Matrices, Second Edition, Academic Press, 1991. 

[31] K A Penson and J-M Sixdeniers, Integral representations of Catalan and related numbers, J. of 

Integer Sequences 4 (2001), Art. 01.2.5 (electronic). 

[32] K A Penson and A I Solomon, Coherent States from combinatorial sequences, preprint 

arXiv:quant-ph/011115ll 2001. 

[33] T Scharf, J-Y Thibon and B G Wybourne, Powers of the Vandermonde determinant and the 

quantum Hall effect, J. Phys. A: Math. Gen. 27 (1994), 4211-4219. 

[34] A Selberg, Bemerkninger om et multiplet integral, Norsk Matematisk Tidsskrift, 26 (1944), 71-78. 

[35] N P Sokolov, Spatial matrices and their applications (in Russian), Gosudarstv. Izdat. Fiz.-Mat. 

Lit., Moscow, 1960. 

[36] N P Sokolov, Introduction to the theory of multidimensional matrices (in Russian), Nukova 





Hankel hyperdeterminants 


32 


Dumka, Kiev, 1972. 

[37] C Radoux, Calcul effectif de certains determinants de Hankel, Bull. Soc. Math. Belg. Ser. B 31 

(1979), 49-55. 

[38] U Tamm, Some aspects of Hankel matrices in coding theory and combinatorics, The Electronic 

Journal of Combinatorics 8 (2001). 

[39] R Vein and P Dale, Determinants and their applications in mathematical physics. Springer, 1999. 

[40] L Vretare, Formulas for elementary spherical functions and generalized Jacobi polynomials, SIAM 

J. Math. Anal. 15 (1984), 805-833. 

[41] N Ja Vilenkin and A U Klimyk, Representation of Lie groups and special functions. Recent 

advances. Mathematics and its Applications, 316. Kluwer Academic Publishers Group, 
Dordrecht, 1995. xvi+497 pp. 



